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Abstract. The spatial distribution of stars in a disk of 
a galaxy can be described by a radial scale length and a 
vertical scale height. The ratio of these two scale parame- 
ters contains information on the axis ratio of the velocity 
ellipsoid, i.e. the ratio of the vertical to radial stellar ve- 
locity dispersions of the stars, at least at some fiducial 
distance from the center. The radial velocity dispersion 
correlates well with the amplitude of the rotation curve 
and with the disk integrated magnitude, as was found 
by Bottema (1993). These relations can be understood 
as the result of the stellar disk being (marginally) stable 
against local instabilities at all length scales. This is ex- 
pressed by Toomre's well-known criterion, which relates 
the sheer in the rotation to a minimum value that the 
radial stellar velocity dispersion should have for stability 
for a given surface density. Via the TuUy-Fisher (1977) 
relation, the velocity dispersion then becomes related to 
the integrated magnitude and hence to the scale length. 
The vertical velocity dispersion relates directly to the scale 
height through hydrostatic equilibrium. It can be shown 
that the ratio of the two length scales relates to the axis 
ratio of the velocity ellipsoid only through the Toomre pa- 
rameter Q and in particular does not require a choice of 
the mass-to-light ratio or a distance scale. 
We have applied this to the statistically complete sam- 
ple of edge-on galaxies, for which de Grijs (1997) has per- 
formed surface photometry and has determined the length 
scales in the stellar light distribution. 

Key words: Galaxies: kinematics and dynamics - Galax- 
ies: photometry - Galaxies: general 



1. Introduction 

It has been known for many decades that the distribution 
of the stellar velocities in the solar neighbourhood is far 
from isotropic. A longstanding problem in stellar dynam- 
ics (or in recent times more appropriately called galactic 
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dynamics) has been the question of the shape and orien- 
tation of the velocity ellipsoid (i.e. the three-dimensional 
distribution of velocities) of stars in disks of spiral galax- 
ies. The local ellipsoid has generated debate and research 
for at least a century. There still is no consensus on the 
question of the orientation of the longest axis of the veloc- 
ity ellipsoid (the "tilt" away from parallel to the plane) at 
small and moderate distances from the symmetry plane of 
the Galaxy (but see Cuddeford & Amendt 1991a, b), which 
is of vital importance for attempts to estimate the local 
surface density of the Galactic plane from vertical dynam- 
ics of stars. The radial versus tangential dispersion ratio is 
reasonably well understood as a result of the local sheer, 
which through the Oort constants governs the shape of the 
epicyclic stellar orbits (but see Cuddeford & Binney 1994 
and Kuijken & Tremaine 1994 for higher order effects as 
a result of deviations from circular symmetry) . 

There are two general classes of models for the origin of 
the velocity dispersions of stars in galactic disks. The first, 
going back to Spitzer & Schwarzschild (1951), is scattering 
by irregularities in the gravitational field, later identified 
with the eflFects of Giant Molecular Clouds (GMCs). The 
second class of models can be traced back to the work 
of Barbanis & Woltjer (1967), who suggested transient 
spiral waves as the scattering agent; this model has been 
extended by Carlberg & Sellwood (1985). Recently, the 
possiblity of infall of satellite galaxies has been recognized 
as a third option (e.g. Velazquez & White, 1999). 

In the solar neighbourhood the ratio of the radial and 
vertical velocity dispersion of the stars CTz/ctr is usually 
taken as roughly 0.5 to 0.6 (Wielen 1977; see also Gomez 
et al. 1990), although values on the order of 0.7 are also 
found in the literature (WooUey et al. 1977; Meusinger et 
al. 1991). The value of this ratio can be used to test predic- 
tions for the secular evolution in disks and perhaps distin- 
guish between the general classes of models. Lacey (1984) 
and Villumsen (1985) have concluded that the Spitzer- 
Schwarzschild mechanism is not in agreement with ob- 
servations: the predicted time dependence of the velocity 
dispersion of a group of stars as a function of age disagrees 
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with the observed age - velocity dispersion relation (see 
also Wielen 1977), while it would not be possible for the 
axis ratio of the velocity ellipsoid CTz/ctr to be less than 
about 0.7 (but see Ida et al. 1993) 

Jenkins & Binney (1990) argued that it is likely that 
the dynamical evolution in the directions in the plane and 
that perpendicular to it could have proceeded with both 
mechanisms contributing, but in different manners. Scat- 
tering by GMCs would then be responsible for the vertical 
velocity dispersion, while scattering from spiral irregular- 
ities would produce the velocity dispersions in the plane. 
The latter would be the prime source of the secular evo- 
lution with the scattering by molecular clouds being a 
mechanism by which some of the energy in random mo- 
tions in the plane is converted into vertical random mo- 
tions, hence determining the thickness of galactic disks. 
The effects of a possible slow, but significant accretion of 
gas onto the disks over their lifetime has been studied by 
Jenkins (1992), who pointed out strong effects on the time 
dependence of the vertical velocity dispersions, in partic- 
ular giving rise to enhanced velocities for the old stars. 

The only other galaxy in which a direct measurement 
of the velocity ellipsoid has been reported, is NGC488 
(Gerssen et al. 1997). NGC488 is a moderately inclined 
galaxy, which enables these authors to solve for the disper- 
sions from a comparison of measurements along the major 
and minor axes. NGC488 is a giant Sb galaxy with a pho- 
tometric scale length of about 6 kpc (in the i?-band) and 
an amplitude of the rotation curve of about 330 km s^^. 
The axis ratio (Tz/ctr is 0.70 ± 0.19; this value, which is 
probably larger than that for the Galaxy, suggests that the 
spiral irregularities mechanism should be relatively less 
important, in agreement with the optical morphology. 

The light distribution in galactic disks has in the ra- 
dial direction an exponential behaviour (Freeman 1970), 
characterised by a scale length h. In the vertical direc- 
tion -at least away from the central layer of young stars 
and dust that is obvious in edge-on galaxies- it turns out 
that the light distribution can also be characterised by an 
exponential scale height /iz, which is independent of galac- 
tocentric distance (van der Kruit & Searle 1981, but see de 
Grijs & Peletier 1997). It then is usually assumed that the 
three-dimensional light distribution traces the distribution 
of mass; this seems justifiable since the light measured is 
that of the old disk population, which contains most of the 
stellar disk mass and dominates the light away from the 
plane. On general grounds, these two typical length scales 
are expected to be independent, the radial one resulting 
from the distribution of angular momentum in the pro- 
togalaxy (e.g. van der Kruit 1987; Dalcanton et al. 1997) 
or that resulting from the merging of pre-galactic units 
in the galaxy's early stages, while the length scale in the 
z-direction would result from the subsequent, and much 
slower, disk heating and the consequent thickening of the 
disk. It is not a priori clear, therefore, that the two scale 
lengths should correlate. Yet, they do bear a relation to 



the ratio of the velocity dispersions of the stars in the 
old disk population. The vertical one follows directly from 
hydrostatic equilibrium. In the radial direction it is some- 
what indirect; a relation between the radial scale length 
and the corresponding velocity dispersion comes about 
through conditions of local stability (e.g. Bottema 1993). 
In a recent study, de Grijs (1997, 1998; see also de 
Grijs & van der Kruit 1996) has determined the two scale 
parameters in a statistically complete sample of edge-on 
galaxies and found the ratio of the two (h/h^) to increase 
with later morphological type. In this paper we will exam- 
ine this dataset in detail in order to investigate whether 
it can be used to derive information on the axis ratios of 
the velocity ellipsoid and help make progress in resolving 
the general issues described above. 

2. Background 

In an extensive study of stellar kinematics of spiral galax- 
ies, Bottema (1993) presented measurements of the stellar 
velocity dispersions in the disks of twelve spiral galax- 
ies. This first reasonably sized sample represented a fair 
range of morphological types and luminosities, although 
it was not a complete sample in a statistical sense. In 
each galaxy he determined a fiducial value for the veloc- 
ity dispersion, namely at one photometric (i?-band) scale 
length. He then found that this fiducial velocity dispersion 
correlated well with the absolute disk luminosity as well 
as with the maximum rotation velocity of the galaxy. His 
sample contained both highly inclined galaxies (where the 
velocities in the plane are in the line of sight) and close to 
face-on systems (where one measures the vertical veloc- 
ity dispersion); when he forced the relations for the two 
classes of galaxies to coincide he found that a similar ratio 
between radial and vertical velocity dispersion as applica- 
ble for the solar neighbourhood was needed. 

Bottema's empirical relation for velocity dispersion 
versus rotation velocity is 



an,h = 0.29 T^ot, 



(1) 



whereas for velocity dispersion versus disk luminosity it 
reads (in the form of absolute magnitude) 



crR,h(km s"^) = -17 X Mb - 279. 



(2) 



These relations can -for any galaxy for which the pho- 
tometry is available or for which the rotation curve is 
known- be used to estimate the radial stellar velocity 
dispersion of the old disk stars at one photometric scale 
length from the center. Doing this for the de Grijs sam- 
ple of edge-on galaxies and estimating the vertical velocity 
dispersion from the vertical scale height, one can in princi- 
ple determine the axis ratio of the velocity ellipsoid for this 
entire sample. It would appear that this is a rather uncer- 
tain procedure, since one will have to assume a mass-to- 
light ratio [M/L) in order to calculate the vertical veloc- 
ity dispersion from the photometric parameters. We will 
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show, however, that M/L does not enter explicitely in the 
formula for the ratio of velocity dispersions. 

We will list our assumptions: 
• The surface density of the disk has an exponential form 
as a function of galactocentric distance: 



I](i?) = S(0) e' 



-R/h 



(3) 



• The vertical distribution of density can be approximated 
by that of the isothermal sheet (van der Kruit & Searle 
1981), but we will use instead the subsequently suggested 
modification (van der Kruit 1988) 



p{R, z) = p{R, 0) sech(z//iz) 



(4) 



A detailed investigation of the sample (de Grijs et al. 1997) 
shows indeed that the vertical light profiles are much closer 
to exponential than to the isothermal solution, although 
the mass density distribution most likely is less peaked 
than that of the light, since young populations with low 
velocity dispersions add significantly to the luminosity but 
little to the mass. Then the vertical velocity dispersion a^ 
can be calculated from 



cr^ = 1.705l7rGI](i?)/iz 



(5) 



The usual parameter Zo used in the notation for the 
isothermal disk (and in de Grijs 1998) is Zo = 2/iz. It is 
important to note that this formula assumes that the old 
stellar disk is self-gravitating. Although this can be made 
acceptable for galaxies like our own at positions of a few 
radial scale lengths from the center (see van der Kruit & 
Searle 1981), it is improbable in late-type galaxies, which 
have significant amounts of gas in the disks, and we will 
need to allow for this. 

• The mass-to-light ratio M/L is constant as a function 
of radius. Support for this comes from the observation by 
van der Kruit & Freeman (1986) and Bottema (1993) that 
the vertical velocity dispersion in face-on spiral galaxies 
falls off with a scale length about twice that of the surface 
brightness (but note that Gerssen et al. 1997, could not 
confirm this for NGC488), combined with the observed 
constant thickness of disks with galactocentric radius. 

• We are not making any assumptions on the functional 
form of the dependence of the radial velocity dispersion 
or the axis ratio of the velocity ellipsoid. The observed 
radial stellar velocity dispersions in Bottema's sample are 
consistent with a drop-off ex exp{—R/2h), in which case 
this axis ratio would be constant with galactocentric dis- 
tance. However, over the range considered the data can be 
fitted also with a radial dependence for the radial veloc- 
ity dispersion in which the parameter Q for local stability 
against axisymmetric modes (Toomre 1964) is constant 
with radius (ex R exp(— i?//i); see van der Kruit & Free- 
man 1986). The definition of Toomre's (1964) parameter 
Q for local stability against axisymmetric modes is 



3.36GE' 



(6) 



Disks are stabilised at small scales through the Jeans 
criterion by random motions (up to the radius of the 
Jeans mass) and for larger scales by differential rotation. 
Toomre's condition states that the minimum scale for sta- 
bility by differential rotation should be no larger than the 
Jeans radius. 

• We assume that spiral galaxies have flat rotation curves 
with an amplitude T^ot over all but their very central 
extent. This assumption implies that we may write the 
epicyclic frequency k as 



:Vr, 



H 

where A and B are the Oort constants. 



(7) 



First we will look into the background of the Bottema 
relations (1) and (2) (see also van der Kruit 1990; Bottema 
1993, 1997). 

Evaluating Toomre's Q aX R = Ih and using the ex- 
pression for the epicyclic frequency above, we find 



3.36G^S(/i)/i 



(8) 



Using S(0) = {M/L)iio and the total disk luminosity from 
id — 27r/io/i^ we get 



1.68G^ fM\ M^iy 



(9) 



Neither in the sample of galaxies that Bottema used to 
define his relations, nor in our sample of edge-on systems 
do we have galaxies with unusually low surface brightness. 
It seems therefore justified to assume that for the galaxies 
considered we have a reasonably constant central surface 
brightness (Freeman 1970; van der Kruit 1987) 



Mo 



21.6 pb = 142 Lq pc^ 



(10) 



where /is stands for 5-magnitudes arcsec"^. So, if ^o, Q 
and [M/L] are constant between galaxies, we see that 
the fiducial velocity dispersion depends only on the disk 
luminosity and the rotation velocity. 

Bottema's relation (1) can then be reconciled with Eq. 
(9), if we have 



Ld ex V^, 



(11) 



This is approximately the TuUy-Fisher relation (TuUy & 
Fisher 1977); not precisely, since we use the disk luminos- 
ity and not that of the galaxy as a whole (however, for 
late-type galaxies this would be a minor difference)F]. 

Eqs. (1) and (2) together would imply an exponential TuUy- 
Fisher relation rather than a power law. The problem is that 
Bottema chose to fit a linear relation to his data, which is 
completely justified in view of his error bars. But eqs. (11) and 
(9) would imply that he should have fitted a curve in which the 
velocity dispersion is proportional to LJ (see van der Kruit 
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So, we see that Bottema's relation (1) follows directly 
from Toomre's stability criterion in exponential disks with 
flat rotation curves as long as Eq. (11) holds. The pro- 
portionality constant in Eq. (11) can be fixed using the 
parameters for the Milky Way Galaxy and for NGC 891 
as given in van der Kruit (1990). These two galaxies have 
Ld - 1.9 X 10^° Lq and V^ot ^ 220 km s~i|. This gives 



Ld (Lq) = 8.11 K*ot (km s-i). 
and using also Eq. (10) we get 

aR,h = 5.08 X 10-2Q fM.) v,,f 



(12) 



(13) 



^From this we find with Bottema's relation (1), that 
Q{M/L)b ~ 5.7. In a somewhat different, but compara- 
ble manner, Bottema (1993) has also concluded that this 
product is of order 5. 

We now turn to the vertical velocity dispersion. Evalu- 
ating the equation for hydrostatic equilibrium (5) at galac- 
tocentric distance R = Ih we find 

a,,h = {5.36 GJ:{h)K}^/^ , (14) 

and can thus calculate the vertical velocity dispersion from 



Cz.h 



^-fcff] Mo/^z 



1/2 



(15) 



Finally we examine the ratio of the two velocity dis- 
persions. If we eliminate S(/i) between Eqs. (8) and (14) 
we obtain 



aR,h = 0.444 Q 



^rot hz 



and with Eq. (1) 






7.77 /iz 



(16) 



(17) 



Note that due to the elimination of the surface density 
also the mass-to-light ratio has dropped out of this equa- 
tion and the result is independent of any assumption on 
M/L. Eq. (17) translates the ratio of the two length scales 
to that of the corresponding velocity dispersions and the 
underlying physics can be summarized as follows. In the 
vertical direction the length scale and the velocity disper- 
sion relate through dynamical equilibrium. In the radial 



1990, p. 199). Performing such a fit to Bottema's data gives 
a curve that is only marginally different from a straight line 
over his range of absolute magnitudes (only a few km s~^). So, 
Eq. (2) should only be seen as an empirical fit of the data to 
a straight line, although these data are equally consistent with 
the dependence following from Eqs. (9) and (11). 

^ The distance for NGC 891 -as are all distances in this 
paper- is based on a Hubble constant of 75 km s~^ Mpc~^. 



direction the velocity dispersion is related to the epicyclic 
frequency through the local stability condition, which is 
proportional to the rotation velocity. The "TuUy-Fisher 
relation" then relates this to the integrated magnitude 
and hence to the size and length scale of the disk. 

One should be careful in the use of Eq. (17), since 
in practice photometric scale lengths are wavelength de- 
pendent and its derivation -and therefore the numerical 
constant- is valid only at one exponential surface density 
scale length. The purpose of presenting it here is only to 
show that, if the two velocity dispersions are derived in a 
consistent manner from Eqs. (8) and (14) -or alternatively 
Eqs. (9) and (15)-, the assumption used for M/L drops 
out in the resulting ratio. 

3. Application to the de Grijs sample 

The sample of edge-on disk galaxies of de Grijs (1998) 
contains 46 systems for which the structural parameters 
of the disks have been determined (including a bulge/disk 
separation in the analysis). From this sample we take those 
for which rotation velocities have been derived in a uni- 
form manner (Mathewson et al. 1992; data collected in de 
Grijs 1998, Table 4) as well as those for which the Galac- 
tic foreground extinction in the 5-band is less than 0.25 
magnitudes. For this remaining sample of 36 galaxies we 
perform the following calculations: 

• From the total magnitudes in de Grijs (1998, Table 6) 
we obtain the integrated B- and /-magnitudes of the disk. 

• Using the radial scale length as measured in the /-band 
we calculate the central (face-on) surface brightness of the 
disk from its /-band integrated luminosity. So, we do not 
use Eq. (10) for a constant central surface brightness. 

• Then we use one or both of the two Bottema relations (1) 
and (2) to estimate the radial velocity dispersion at one 
photometric scale length (by definition in the _S-band). 
Where we can do this with both relations, the ratio be- 
tween the two estimates is 1.11 ± 0.19. This is not trivial, 
as the rotation velocities and disk luminosities are deter- 
mined completely independently (and by different work- 
ers) and only the one using the absolute magnitude needs 
an assumption for the distance scale. 

• Then using Eq. (15), we estimate the vertical velocity 
dispersion at one photometric scale length in the /-band. 
For this we need a value for the mass-to-light ratio and 
we will discuss this first. 

We found, that through Eq. (13) Bottema's relation 
(1) provides a value for Q{AI/L) of about 5.7. So we make 
a choice for Q rather than for M/L. It has become cus- 
tomary to assume values of Q of order 2, mainly based on 
the numerical simulations of Sellwood & Carlberg (1984), 
who find their disks to settle with Q ~ 1.7 at all radii. In 
principle we can use the observed properties of the Galaxy 
to fix Q from Eq. (17). We have (crz/crR)^ ~ 0.5 (in the 
solar neighbourhood, but assume for the sake of the argu- 
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ment also a.t R = Ih) and h^/h ^ 0.1 (see Sackett 1997 
for a recent review), so that indeed Q ~ 1.7. We will make 
the general assumption that Q = 2, in agreement with the 
considerations above; then {M/L)b = 2.8. 

The rotation velocity version of Bottema's empirical 
relations (Eq. (1)) can provide further support for the 
choice of Q, along the lines of the discussion in van der 
Kruit & Freeman (1986). In the first place we recall the 
condition for the prevention of swing amplification in disks 
(Toomre 1981), as reformulated by Sellwood (1983) 



X 



Rn^ 



27rmGI] 



>3, 



(18) 



where m is the number of spiral arms. For a flat rotation 
curve this can be rewritten as 



QVro 

CTR 



> 3.97TO. 



(19) 



With Eq. (1) this becomes Q > 1.15to. Considering that 
the coefficient in Eq. (I) has an uncertainty of order 15%, 
this tells us that we have to assume Q at least of order 2 
to prevent strong barlike (m=2) disturbances in the disk. 
A similar argument can be made using the global sta- 
bility criterion of Efstathiou et al. (1982). This criterion 
states that for a galaxy with a flat rotation curve and an 
exponential disk, global stability requires a dark halo and 



Y = V,, 






1/2 



1.1. 



(20) 



Here M^^^k is the total mass of the disk. This can be rewrit- 
ten as 



Y = 0.615 



QRV,c 



han 



1/2 



exp 



R 

2h 



1.1, 



(21) 



and, when evaluated at i? = l/i, yields with Eq. (1) 
0.69\/Q ^ 1.1, and therefore also imphes that Q should 
be at least about 2. Efstathiou et al. have also come to 
this conclusion for our Galaxy, with the use of local pa- 
rameters for the solar neighbourhood. 

Having adopted a value for Q and through this a value 
for {M/L)b, we will have to convert it to {M/L)j. For this 
we need a. {B — I) colour for the disks. From the fits of de 
Grijs (1998) we find that the total disk magnitudes show a 
rather large variation in colour; for the sample used here 
{B — I) has a mean value of 1.9, but the r.m.s. scatter 
is 0.8 magnitudes. In his discussion, de Grijs (1998) sus- 
pects a systematic effect of the internal dust in the disks 
(particularly on the 5-magnitudes, which is another rea- 
son for us to use the Vrot -version of the Bottema relation 
(Eq. (1)) in our derivation in the previous section). In- 
stead we turn to the discussion of de Jong (1996b), who 
compares his surface photometry of less inclined spirals to 
star formation models. From his Table 3, we infer that for 
single burst models with solar metallicity and ages of 12 



Gyr {M/L)b = 2{M/L)i. So we will use an (i\//i)/ of 
1.4. 

There is a further refinement required. In order to take 
into account the fact that in late-type galaxies the gas con- 
tributes significantly to the gravitational force, we have to 
correct for a galaxy's gas content as a function of Hubble 
type. In the following, we will discuss the observational 
data regarding the Hi and the H2 separately. 

For 25 of de Grijs' sample galaxies Hi observations are 
available, so that we can estimate the gas-to-total disk 
mass. For this we apply a correction of a factor 4/3 to the 
Hi in order to take account of helium and use de Grijs' 
(1998) /-band photometry and our adopted M/Lb ratio 
of 2.8 (see below) to estimate the total disk mass. As a 
function of Hubble type we then find 



Type gas-to-total disk mass n 



Sb 


0.31 ± 0.17 


3 


Sbc 


0.36 ±0.19 


5 


Sc 


0.53 ± 0.09 


5 


Scd 


0.49 ± 0.16 


9 


Sd 


0.52 ± 0.10 


3 



We find no dependence on rotation velocity: 


l^-ot (km s~^) gas-to-total disk mass 


n 


80 - 130 0.51 ± 0.09 
130 - 180 0.45 ± 0.17 
180 - 230 0.51 ± 0.09 


10 

8 

7 



So, the Hi mass is about half the stellar mass in disks of 
Sb's and about similar to that in Sc's and Sd's. But there 
is no dependence on rotation velocity. But this is not what 
we need; we should use surface densities rather than disk 
masses. Now, the Hi is usually more extended than the 
stars and has a shallower radial profile. So the ratios in 
the tables above are definite upper limits. In order to take 
into account the effect that in late-type systems the gas 
contributes significantly to the gravitational force we have 
"added" for types Scd and Sd a similar amount of gas as 
in stars and half of that for Sc's. 

The distribution of H2 in spiral galaxies is a more com- 
plex matter; it is often centrally peaked, although some 
Sb galaxies exhibit central holes (for a recent review, see 
Kenney 1997). The molecular fraction of the gas appears 
to be lower in low-mass and late-type galaxies, assuming 
that the conversion factor from CO to molecular hydrogen 
is universaln. Since our sample galaxies are generally low- 
mass, later-type systems, we believe that the corrections 
for molecular gas are small, and therefore contribute little 
to the correction for the presence of gas. 

We added (a) the galaxies from van der Kruit & Searle 
(1982) to the sample, (b) our Galaxy using the Lewis & 



^ There is even doubt concerning the constancy of this con- 
version factor within our Galaxy (Sodroski et al. 1995). 
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Fig. 1. The scale lengths (a) and scale heights (b) of the galaxies in our sample as a function of their rotation velocities. 
The r.m.s. errors are of order 10 km s~^ in the rotation velocities (Mathewson et al. 1992), 5% in the radial scale 
length and 9% in the vertical scale height (de Grijs 1998). 



Freeman (1989) velocity dispersion and the structural pa- 
rameters in van der Kruit (1990), and (c) the observational 
results for NGC488 from Gerssen et al. (1997). We leave 
the few early type (SO and Sa) galaxies out of the dis- 
cussion, because the component seperation in the surface 
brightness distributions is troublesome and some of our 
assumptions (in particular the self-gravitating nature of 
the disks) are probably seriously wrong. 

In order to be able to trace the origin of our results, 
we first show in Fig. 1 the radial and vertical scale lengths 
of the sample as a function of the rotation velocity. Both 
increase with K-ot) which would be expected intuitively. 
The main result is presented in Figs. 2 and 3. From Fig. 
2 we see that the vertical velocity dispersions, that have 
been derived from hydrostatic equilibrium, increase with 
the rotation speed (the radial velocity dispersions do the 
same automatically as a result of the use of the Bottema 
relations). For the slowest rotation speeds the predicted 
vertical velocity dispersion is on the order of 10-20 km s^^, 
which is close to that observed in the neutral hydrogen in 
face-on galaxies (van der Kruit & Shostak 1984). 

The distribution of the axis ratio of the velocity ellip- 
soid with morphological type is as follows: 



Type 


0'z.h/o'R„h 


n 


Sb 


0.71 ± 0.14 


11 


Sbc 


0.69 ± 0.16 


7 


Sc 


0.49 ± 0.17 


6 


Scd 


0.70 ± 0.20 


11 


Sd 


0.63 ± 0.22 


5 



Not much of a trend is seen here. It is in order to comment 
here briefly on the effects of our corrections for the gas to 
obtain vertical velocity dispersions. From our discussion 
above we conclude that there would be no systematic ef- 
fect introduced as a function of rotation velocity. Further- 
more, taking away our correction altogether reduces the 
values for the average axis ratio in the table just given to 
about 0.55 for Scd and Sd galaxies. Even in this unreal- 
istic case of not allowing for the presence of the gas, we 
believe the trend to be hardly significant in view of the 
uncertainties. Since some correction for the gas mass as 
a function of morphological type must be made, we can- 
not claim that we find any evidence for a change in the 
velocity anisotropy with Hubble type. 

Fig. 3 shows the axis ratio of the velocity ellipsoid of 
all the galaxies versus their rotation velocity. In view of 
the fact that the two dispersions that go into this ratio are 
determined from different observational data (cR^h from 
integral properties such as total luminosity and amplitude 
of the rotation curve; az^h from photometric scale parame- 
ters and surface brightness) and that we have made rather 
simplifying assumptions, the scatter is remarkably small. 
No systematic trends are visible (and would probably not 
be significant!) in the data. The points closest to unity 
in the dispersion ratio generally have low rotation veloc- 
ities and inferred velocity dispersions. One of these (K-ot 
= 95 km s-i, crz,h/CTR,h = 0.75) is NGC5023. Bottema et 
al. (1986) have shown that the stars and the gas in this 
galaxy are effectively coexistent; the radial and vertical 
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Fig. 2. The calculated vertical velocity dispersions at one exponential scale length from the center as a function of the 
rotation velocities 



distributions are very similar. This would imply that the 
velocity dispersions of the gas and the stars are the same. 
The vertical velocity dispersion found here is about 20 
km s~^, which is significantly higher than that observed 
in larger spirals. It would be of interest to measure the Hi 
velocity dispersion in this galaxy. Since the HI would be 
expected to have an isotropic velocity distribution from 
collisions between clouds, the vertical dispersion should 
be equal to that in the line of sight in edge-on galaxies. 



4. Discussion 

In this section we will critically discuss the uncertainties 
in our approach. 

• The linearity of the magnitude version of the Bot- 
tema relation. Wc have discussed above that the power- 
law nature of the TuUy-Fisher relation (Eq. (11)) would 
imply a nonlinear form of the magnitude version of Bot- 
tema's relation (Eq. (2)). We have used it as an empirical 
relation to help (together with Eq. (1)) to estimate the 
radial velocity dispersions from the observed photometry. 
One may argue that it is internally consistent to use in- 
stead of Eq. (2) a fit of the form 



rl/4 



(22) 



This has only a noticable effect on galaxies with faint abso- 
lute disk magnitudes. We have repeated our analysis using 
such a fit and find no change in our results. To be more 
definite we repeat the table of the average axis ratio as a 
function of morphological type that we then obtain. 



Type 


0-z,h/0'R,h 


n 


Sb 


0.67 ± 0.17 


11 


Sbc 


0.64 ±0.13 


7 


Sc 


0.47 ±0.17 


6 


Scd 


0.60 ± 0.11 


11 


Sd 


0.60 ± 0.20 


5 



• Choice of Q. We have adopted a value for Q of 2.0, 
and this enters directly in our results in the calculation of 
the vertical velocity dispersion through the value of M/L 
that follows from this choice. Had we adopted a value of 
1.0 for Q, M/L would have been a factor 2 higher and the 
value for a^^h a factor ^/2 -see Eq. (15)-. Bottema (1993, 
his Fig. 11) has shown that his observations of the stellar 
kinematics do not show any evidence for systematic vari- 
ations in Q among galaxies. 

We have used Eq. (19) to argue that Q is likely of order 
2 in order to prevent strong barlike (m=2) disturbances. 
It does not necessarily follow from this that galaxies with 
more spiral arms should have higher values of Q or that Q 
should be significantly lower in galaxies with very strong 
two-armed structure. Spiral structure may arise in a vari- 
ety of ways; we only argue that disks do not have grossly 
distorted m=2 shapes and that therefore swing amplifica- 
tion is apparently not operating. 

We have assumed that the stellar disk is self-gravitating 
and ignored the influence of the gas in the evaluation of 
Q. At one scale length this is probably justified (it is only 
a small effect in solar neighbourhood), even for late- type 
disks. 
• Non- exponential nature of the disks. Often the disk in 
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Fig. 3. The ratio of the calculated vertical to radial velocity dispersions at one exponential scale length from the center 
as a function of the rotation velocities 



actual galaxies can be fitted to an exponential only over 
a limited radial extent. In that case our description is un- 
likely to hold. However, in our sample the fits can be made 
reasonably well at one scale length from the center and we 
believe this not to be a problem. 

• Vertical structure of the disks. Our results depend on the 
adoption of a particular form of the vertical mass distribu- 
tion (namely the sech(z)-form) of Eq. (4). This enters our 
results through the value of the numerical constant 1.7051 
in Eq. (5), and in Eq. (15) it enters into the value for a^^h 
as its square root. Had we assumed the isothermal dis- 
tribution, then the constant would have been 2.0, while 
it would have been 1.5 for the exponential distribution. 
This would have given us values for a^^h which are only 6 
to 8% higher or lower. As we have shown in de Grijs et 
al. (1997), for our sample galaxies the vertical luminosity 
distributions in these disk-dominated galaxies are slightly 
rounder than or consistent with the exponential model. 
However, the vertical mass distribution is probably less 
sharply peaked, and thus expected to be more closely ap- 
proximated by the sech(2) model. 

• Non-constancy of central (face-on) surface brightness. 
We have assumed in section 2 that for all galaxies the cen- 
tral surface brightness is constant. This is certainly unjus- 
tified for so-called "low surface brightness galaxies" ; how- 
ever, our galaxies have brighter surface brightnesses than 
galaxies that are usually considered to be of this class. But 
even for galaxies as in our sample it remains true that the 
(face-on) central surface brightness is in general somewhat 
lower for smaller systems (van der Kruit 1987; de Jong 



1996a). From de Jong's bivariate distribution functions, it 
can be seen that late-type, low absolute magnitude spi- 
rals may have a central surface brightness (in B) that is 
up to 1.0 magnitude fainter. Note, however, that for our 
derived velocity dispersons we have used the actually ob- 
served surface brightness for each galaxy. 

• Effects of colour variations on the mass-to-light ratio. 
There is a fairly large variation in the colours of the disks 
in the sample. De Grijs (1998) has argued that this is the 
result of internal dust extinction. However, we have used 
the /-band data; de Grijs' Fig. 11 shows that the varia- 
tion is much less in (/ — K) than in colours involving the 
-B-band. 

The colours of the disks do not correlate with morpholog- 
ical type (de Grijs 1998); although such a correlation has 
been seen in more face-on galaxies (de Jong 1996b), we 
believe that the fitting procedure has ignored most of the 
young population and dust absorption near the plane, and 
that contributions to the M/L scatter as a result of young 
populations has mostly been avoided; and that the /-band 
scale lengths determined away from the galactic planes are 
fairly representative of the stellar mass distributions (de 
Grijs 1998). 

The colour of the fitted disks in our sample has (/ — K) 
in the range ^ 2 to 4. The latter is red, even for an old 
population and may well be caused by excessive internal 
extinction, but we see no strong evidence for a substantial 
systematic correction in our velocity dispersions from this. 

• Effects of metallicity on the mass-to-light ratio. De Jong 
(1996b) has drawn attention to the non- negligible effects 
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of metallicity on the mass-to-light ratio. From his compi- 
lation of models, in particular his W94 (Worthey 1994) 
models with ages of 12 Gyr, we estimate that the effect 
in the /-band amounts to 10 to 20% over the range of 
relevant metallicities. The effect on the derived velocity 
dispersions is the square root of this. 

• Effects of radial colour variations. Since we use an em- 
pirical relation to derive the radial velocity dispersion at 
one -S-scale length from the center, we have to consider 
the effect of using the /-band. We have used the latter as 
the proper scale length to use for the mass density distri- 
bution. The correct one to use here would be the one mea- 
sured in the /f-band, which is 1.15 ± 0.19 times smaller 
for this sample. The scale length in the S-band is 1.64 ± 
0.41 times longer than in the /f-band (values quoted here 
from de Grijs 1998). We deduce from this that we may 
have underestimated the scale length to use by a factor 
1.43 and systematically overestimated the vertical veloc- 
ity dispersion by about 20%. 

The effect of radial metallicity variations on the scale 
length is probably very small; these gradients in the older 
stellar populations are in any case expected to be signifi- 
cantly less than in the interstellair medium. This is so, be- 
cause in models for galactic chemical evolution the mean 
stellar metallicity of the stars approximates the (effective) 
yield, while that in the gas grows to much large values in 
most models (van dcr Kruit 1990, p. 322). 

• Non-flatness of rotation curves. This may be an effect 
for small, late-type galaxies that have slowly rising rota- 
tion curves. It enters however in our analysis only in the 
derivation of the Bottema relations and this holds empir- 
ically to rather small rotation velocities (^ 100 km s^^). 

• The slope of the Tully-Fisher relation. Although we use 
a relation between the luminosity of the disk alone and 
the rotation velocity, it remains true that our "slope" of 
Vcoj is steeper than the usually derived slopes of Tully- 
Fisher relations, which would indicate y^^j (e.g. Giovanelli 
et al.'s 1997 "template relation" yields an exponent of 3.07 
± 0.05). If this slope were to be put into Eq. (9), Eq. (1) 
would have CTR.h ex y/Vmt, which is very significantly in dis- 
agreement with Bottema's observations. The same holds 
for the other Bottema relation, Eq. (2). 

The analysis of the present sample (see de Grijs & Peletier 
1999) has resulted in slopes in the Tully-Fisher relation 
of 3.20 ± 0.07 in the /-band and 3.24 ± 0.21 in the K- 
band. On the other hand, Verhcijen (1997), in his exten- 
sive study of about 40 galaxies in the Ursa Major cluster, 
finds a slope of 4.1 ± 0.2 in the /f-band. 

• Effects of the gas on the value of Q. In our calculations 
we have already made crude allowance for the effects of 
the gas on the gravitational field. But in our derivations 
we have not taken into account the effect of the Hi on the 
effective velocity dispersion to be used in the evaluation 
of the Q-paramcter. The effect of the Hi is a decrease of 
the effective velocity dispersion and therefore in Q. This 
means that the assumed value should in reality be de- 



creased on average, but beyond this numerical effect, it 
does not affect our results. 

The effects just discussed can produce errors in the 
estimated velocity dispersions of the order of 10 to 20% 
each. The final result of the dispersion ratios in Fig. 3 may 
therefore be wrong by a few tenths, which is comparable 
to the scatter in that figure. However, we have no cause to 
suspect that we have introduced serious systematic effects 
that would be strong functions of the rotation velocity or 
the morphological type and the lack of correlation of the 
axis ratio with these properties is unlikely to be an artifact 
of our analysis. 

We conclude that it is in principle possible to infer 
information on the axis ratio of the velocity ellipsoid from 
a sample of edge-on galaxies for which both the radial 
scale length as well as the vertical scale height have been 
measured. The result, however shows much scatter, most 
of which is a result of the necessary assumptions. There 
is one significant improvement that can be made and that 
is the direct observation of the stellar velocity dispersion 
in these disks. That this is feasible in practice for edge-on 
systems has been shown by Bottema et al. (1987, 1991). 
The observed velocity profiles can be corrected for the line- 
of-sight effects, giving the tangential velocity dispersion, 
which through the observed shape of the rotation curve 
can be turned into the radial velocity dispersion. Although 
a time consuming programme, we believe that it is worth 
doing for two reasons: (1) It will set both versions of the 
Bottema relation on a firmer footing. (2) The uncertainties 
in the analysis above can likely be significantly diminished 
by direct measurement of the radial velocity dispersion 
rather than having to infer it from the rotation velocity 
or the disk absolute magnitude. 
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